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Abstract 

In this paper we try to answer the main question: what is a quan- 
tum black hole? 
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1 Introduction 



What does urge a researcher to investigate quantum black holes? Honestly, 
his own intrinsic interest. Besides, there are other, more (or less) scientific 
reasons. It is commonly believed that only small black holes can be consid- 
ered as quantum objects. Small, what does it mean? To estimate, we should 
compare the size of the black hole with the corresponding Compton length. 
The gravitational radius Vg of the black hole of mass m equals Vg = 
where G is the Newtonian gravitational constant, and c is the speed of light. 
The Compton length of such particle is A = ^ is the Planck's constant). 

If Tg ~ A, than the so called Planckian mass is mpi = ~ 10"^^^, the 

Planckian length is Ipi = \J^- ~ IQ^^^cm. In what follows we us the units 

in which ^ = c = 1, so nipi = 1/y/G, Ipi = \fG. Black holes of so small mass 
and sizes could be created from large metric fiuctuation in the very early 
Universe (primordial black holes), or during vacuum phase transition. 

The mankind "knows" about black holes since 1784, when some priest J. 
Michell recognized that stars becomes invisible if the correspond escape 
velocity will exceed the speed of light. He called such stars "dark stars" 
Fifteen years later, in 1799, P.S. Laplace derived a relation between the 
radius and the mass of these dark stars 0. 

Let us consider a spherically symmetric gravitating body of mass M and 
of radius Rq. Then the escape velocity Ve can be determined from the balance 
of the kinetic and potential energies of the particle that starts to move away 
from this gravitating body. Using the non-relativistic mechanics and non- 
relativistic Newtonian theory of gravity we obtain, following Laplace [G is 
the Newton constant, and m is the mass of the escaping particle) 

mvl GmM 
^ Ro ' 

(1.1) 

2 _ 2GM 

Equating the escape velocity to that of light c we get for the maximal radius 
of the invisible body 

_ 2GM 

Ro — n • 
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We are used to think of the black holes as of the extremely dense bodies. 
Indeed, the black hole with the mass of the Earth has the radius about 3 cm 
and the density much higher than the nuclear density. But due to the linear 
radius-mass relation the mean density p is proportional to the inverse mass 
squared, 

3c^ 

^ ~ 327rG'3M2 ■ 

And for the huge black hole formed, say, of the cluster of galaxies, the density 
may be less than that of the atmospheric air. Thus the observer falling 
into such a black hole would not recognize that he never would go out of it 
(moreover, the free fall time in this case is about 10,000 years). 

The history of the "real" black holes started in 1916, when K.Schwarzschild 
found the, now famous, solution to the Einstein equation for the spherically 
symmetric vacuum space-time outside gravitating point source 0. He de- 
rived it in the form 

ds" = Fdt" - F-'dr^ - r^dn^ (1.2) 

(1.3) 

F = l-^ 

r 

Here dfl^ is the line element of the unit sphere, r - its radius, and m is 
the total mass (energy) of the system, also known as Schwarzschild mass. 
The line element ( p..2|) has a singularity at the so called Schwarzschild radius 
Tg = 2Gm, from where the light cannot escape to the infinity. It is surprising 
that value of the Schwarzschild radius calculated using General Relativity 
coincides exactly with non-relativistic (Newtonian) radius of the dark stars. 

In classical General Relativity black holes are very special (and, there- 
fore, very interesting) objects. First of all, they are universal in the sense 
that they are described by only few parameters: their mass, angular momen- 
tum and charges. In the process of black hole formation, (i.e., in the process 
of gravitational collapse) all higher momenta and non conserved charges are 
radiating away. This feature is formulated as a following conjecture: "black 
holes has no hairs". Thus, a black hole formation results generically in the 
loss of information about initial states and previous history of collapsing mat- 
ter. The boundary of the black hole, the so-called event horizon, is the null 
hyper-surface that acts as a one-way membrane. The matter can fall inside 
but can not go outside. Because of this the area of black hole horizon can 



3 



not decrease. These two features, the loss of information and nondecrease of 
the horizon area, allowed J.Bekenstein to suggest the analogy between the 
black hole physics and thermodynamics and identify the area of the horizon 
with the entropy (up to some factor). He did this for the simplest, spher- 
ically symmetric neutral (Schwarzschild) black hole which characterized by 
only one parameter, Schwarzschild mass. Later the four laws of thermody- 
namics were derived for a general black hole. 

In thermodynamics the appearance of the entropy is accompanied by the 
temperature. While the nature of the black hole entropy was more or less 
clear, the notion of its temperature remained mysterious until the revolu- 
tionary work by S. Hawking |^. He showed that the black hole temperature 
introduced by J.Bekenstein is the real temperature in the sense that the black 
hole radiates, and this radiation has a Planck's spectrum. The entropy ap- 
peared equal one fourth of the event horizon area divided by Planckian length 
squared. Thus, even large (compared to the Planckian mass and size) black 
holes exhibit quantum features. It should be stressed that such quantum 
effect is global, namely, it emerges as a result of nontrivial boundary condi- 
tions for the wave function of a quantum field theory in curved space-times 
nontrivial causal structure (existence of the event horizon(s)). 

Due to the process of Hawking's evaporation any (even super-large) black 
hole becomes eventually small enough to be considered as a (local) quantum 
object. At this stage the combination of the global and local quantum effects 
may result in unpredictable features. That why it is so exciting to try to 
understand quantum black hole physics. 

2 Classical Black Holes 

Everybody knows what the classical black hole is. In short, black hole is a 
region of a space-time manifold beyond an event horizon. In turn, an event 
horizon is a null surface that separates the region from which null geodesies 
can escape to infinity and that one from which they cannot. The most 
general black hole is characterized by only few parameters, its mass, angular 
momentum, and some conserved charge, similar to the electric one. In general 
the space-time outside a black hole is stationary and axially symmetric, but 
for a non-rotating black hole it is spherically symmetric and static. In what 
follows we are dealing only with spherically symmetric black holes. It is 
important to stress that the notion of the of the event horizon is global, it 
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requires knowledge of both past and future histories. To understand this 
better, let us consider the general structure of spherically symmetric space- 
time manifolds. The metric of such space-time can always (at least locally) 
be written in the form 

^ gQQdf + 2goidtdq + gudq^ - R^dVL^ (2.1) 

Four metric coefficient are functions of some time coordinate t and some 
radial coordinate q, dVL = 9"^ + sin6'^d(f)'^ and R(t,q) is the radius in that 
sense that the sphere area is AirR'^. Making use of two allowed coordinate 
(gauge) transformations it is always possible to transform a two-dimensional 
(t,q) - part of the metric to the conformally fiat metric. Hence any spherically 
symmetric space-time is (locally) described, actually, by only two invariant 
functions. One of them is, evidently, the radius R{t, q). For the other one it 
is convenient to choose the squared normal vector to the surfaces of constant 
radius 0, i 

A = g''^R,iRj (2.2) 

The function A(t, q) brings a nontrivial information about a space-time struc- 
ture. Indeed, in the fiat Mimkowskian space-time A = 1, all the surfaces 
R = const, are time-like and therefore, R can be chosen as spatial coordi- 
nate q = R on the whole manifold. But in the curved space-time A is no 
more constant and can in general be both positive and negative. The region 
with A > is called the i?-region, and the radius can be chosen as a radial 
coordinate q. In the region with A < the surfaces R = const, are space-like 
(the normal vector is time-like), and the radius R can be chosen as a time 
coordinate t. Such regions are called the T-regions. The R- and T-regions 
were introduced by Igor Novikov. But this is not the whole story. It is easy 
to show that we can not get -R = ("dot" means a time derivative) in a 
T-region. Hence it must be either R > (such a region of inevitable expan- 
sion is called T+-region) or i? < (inevitable construction, a T_-region. The 
same holds for i?-regions. They are divided in two classes, those with R' > 
("prime" stands for a spatial derivative) which are called -R+-regions, and 
i?_-regions with R' < 0. These, R- and T-regions are separated by the sur- 
faces A = which are called the apparent horizons. The apparent horizons 
can be null, time-like or space-like. If surface R = lies inside a T_-region, 
then boundary of the latter forms the so-called trapped surface from which 
any light can escape to the infinity. This part of the apparent horizon lies 
either inside a black hole or coincides partly with an event horizon. Clearly, 
the notion of the apparent horizon is local. 
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After all these preparation we consider the following example. Let us 
imagine that we have a trapped surface and, hence a black hole , and spheri- 
cally symmetric matter layer outside it. let such a layer is held at some fixed 
value of radius with the help,say,of rocket engine. If the engine is abruptly 
switched off, the layer will start to collapse and, eventually, will cross the 
apparent horizon. Again a black hole will be formed, but now it will have 
different parameters (higher mass), and a new event horizon has to be de- 
fined. In this case some part of the apparent horizon (that one lying in the 
past light cone of the crossing point ) will appear to be inside the new black 
hole, while the other part (the new one lying in the future of the crossing 
point) will coincide with the new event horizon. In principle, we may have 
many different matter layers, and nobody can foresee, when the astronauts 
will decide to switch off engines. Our simple gedanken experiment demon- 
strates clearly the global nature of the event horizon. 

All the notions introduced above are extremely important in studying 
black holes formed during gravitational collapse. But there exist a rather 
simple and, nevertheless, very important model, called the eternal black hole 
which we study in the next Subsection. 



2.1 Schwarzschild black hole 



All the Schwarzschild metrics ( |1.2|) forms the one-parameter family labeled 
by the value of total mass m, and there are no other spherically symmet- 
ric solutions to the vacuum Einstein equation which asymptotically flat at 
infinity (this statement is known as Birkhoff's theorem). It is clear at once 
that the metric ( |1.2| ) has a singularity at the gravitational radius Vg = 2Gm. 
The nature of this singularity was not understood for a long time. And only 
after works by M.Kruskal |T0[ and I.D.Novikov it become clear that it is 
actually a coordinate singularity, appeared due to the static form of the line 
element ( |1.2| ) and impossibility to synchronize the clocks of static observers 
at spatial infinity (r = const oo) and of the observers (which are inevitable 
non- static) in the region r < 2Gm. 

The manifold, described by the line element (|1.2|) , is not geodesically com- 
plete. By definition, in geodesically complete space-time (which also called 
a maximally analytically extended manifold) any time-like or null geodesies 
should start and end either at the infinity or the singularity. The structure of 
maximally extended spherically symmetric manifold very easy to see and in- 
vestigate using the so-called Carter - Penrose conformal diagrams. On these 
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diagrams all the infinities lie at the finite distance from the center, and every 
point represents a sphere. 

The conformal diagram of geodesically complete Schwarzschild space-time 
is shown in Fig.l. 




Figure 1: 

We see that this manifold has a geometry of the non-traversable worm- 
hole. There arc two asymptotically flat i?-|--regions connected by the Einstein- 
Rosen bridge. If we start from the left spatial infinity at and go in the 
direction to the right infinity at i^, the radii of spheres first decrease in the 
i?_-region, reach minimal value on the Einstein-Rosen bridge and then in- 
crease in the i?+-region. The narrowest part of the bridge called "throat" , it 
lies at the intersection of two null surfaces r = 2Gm, which are the horizon. 
For the sake of simplicity we will often refer to the i?+-region as the 'our 
side', and to the i?_-region as the 'other side' of the Einstein- Rosen bridge. 
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The gravitating sources in the eternal Schwarzschild black hole are assumed 
to be concentrated in the future and past T-regions at the singular space-like 
surfaces of zero radius, R = 0. These singularities are the real ones, because 
(as can be easily shown) the Riemann curvature tensor and, hence the tidal 
forces, become infinite here. The Schwarzschild geometry in the equatorial 
plane at the fixed moment of time can also be visualized a geometry of some 
surface of rotation embedded into three-dimensional fiat space, (see Fig. 2). 
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Figure 2: 

In what follows we will deal sometimes with the electrically charged 
sources. The space-time outside such sources is described by the Reissner- 
Nordstrom line element which differs from the metric (|1 . 2|) by different func- 
tion F, now it becomes F(r) = 1 — + where e is the electric charge. 

2.2 Sources 

Eternal black hole is the nice and transparent model. It appeared very useful 
in investigating of classical black holes. But it really physical because of the 
absence of any dynamical degrees of freedom, namely, there's no dynamics of 
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gravitating matter, neither real mass-shell spherically symmetric gravitons. 
Therefore, having in mind the construction of quantum black hole models, we 
have to include more realistic gravitating source into consideration. Because 
of the nonlinearity of the Einstein equations the inclusion of source make the 
quantization problem very complicated and, in general, unsolvable. There- 
fore, in order to get some definite results we have to choose the most simple 
types of sources. We will consider the spherically symmetric thin dust shell 
which is a direct generalization of a point gravitating mass. 

The theory of thin shells in General Relativity was developed by W. Israel 
|Tl| and applied then by many authors studying problems of gravitational 



collapse and cosmology [|T2, 14]. We don't need yo consider this theory 
in details. In our case, both outside and inside the shell we have either 
vacuum (Schwarzschild) or electro- vacuum (Reissner- Nordstrom) metric, and 
and complete information about the shell dynamics is contained in the single 
equation, namely, in the (Q)-Einstein equation for the shell (this is the so 
called energy constraint). It has the form 

C^m\/p^ + - Gout^ + Ao„f = ^'KGpSl (2.3) 

Here p is the shell's radius as a function of the proper time of the observer 
sitting on the shell, the dot means the proper time derivative, is the 
surface energy density concentrated on the shell, and Aj„ is the value of the 
(introduced above) invariant function just inside and outside the shell. In our 
case it equals A = —F. The step function A take two values ±1 depending 
om whether radii increase in the direction of outward shell's normal vector 
(cr = +1) or they decrease (a = — 1). It is clear that the sign of a coincides 
with the sign of the i?-region where the shell moves, and its sign can be 
changed only in the T-regions. For the sake of simplicity we consider here in 
details only the case when the shell is not electrically charged and there is 
no other gravitating sources inside it. Then 

Ajri, 1, ^ in. ^ 1 



^out = 1 - ^ + ^, m > 



(2.4) 



Our shell is made of dust, that is, of noninteracting (except gravitationally) 
particles. This means that th surface tension is zero, S\ = S\ = 0. Then, 
from the continuity condition, 

^ + ^(5°-5,^) + [ToT=0 (2.5) 
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it follows that 

where M is readily identified as a bare mass of the shell (the sum of the rest 
masses of particles without the gravitational mass defect). 

Finally, we arrive at only one equation which brings all the information 
about the shell's motion, namely 



r-, / , 2Gm GM , 

Vp2 + 1 - ooutslp' + 1 = 2.7 

V p p 

We can solve this equation for the total mass m to get 

m = M^Jp" + 1 - -—- (2.8) 
2p 

and now it is seen that it is nothing more but the energy conservation equa- 
tion, the square root being the famous Lorentz factor written in terms of 
proper time derivatives. 

We shall be interested here in bound motion only. Let us denote by po 
the radius of the shell at the moment of rest, and let po > R+ (that is, we 
are outside the event horizon). Then it follows from Eqn.(p^) and Eqn. ( p^ 
that m < M and 

^ = M-^^ (2.9) 



2po 

aout = Stgn{l 2.10 

V Po 



Differentiating Eqn. (|2.9 | ) with respect to the bare mass M we obtain 

dm _ GM 
dM~ " ^ 



(2.11) 



Let us consider the dependence of the total mass m on the bare mass M 
for some fixed value of the turning radius po. As seen form Fig. 3, the curve 
has two branches, the increasing and the decreasing ones. 
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Figure 3: 
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The total mass m first increase with M, and the parameters of the shell 
in this case are such that 

— > 

(Jout = +1 (2.12) 

i < — < 1 

The world lines of these shells start from the past singularity at -R = 0, pass 
the i?+-region and end at the future singularity at i? = 0. We shall call such 
a situation the black hole case (the shell forms the black hole) . The complete 
Carter-Penrose diagram is shown in Figs.4,5. 

The spatial geometry at the moment of time symmetry is shown in Fig.6. 



The curve m{M) reaches then the maximum and starts to decrease. For 
the decreasing branch we have 

— < 

(Tout = -1 (2.13) 

m ^ 1 
M ^ 2 

Now the shells pass the i?_-region, and we shall call this the wormhole case 
(the shell forms the wormhole on the other side of the Einstein- Rosen bridge). 
The corresponding Carter-Penrose diagrams is shown in Fig. 8, 9. The spatial 
geometry is shown in Fig. 7. It is explicitly seen in the figure that the radii 
decrease outside the shell, reach the minimal value 2Gm at the bifurcation 
point and then they start to increase. The wormhole region can not be 
reached or seen from the i?+-region. In the marginal case m/M — 1/2 and 
the turning point lies exactly at the bifurcation point R — 2Gm. 

It is instructive to consider some limiting cases. The first one is the limit 
M ^ oo. When the bare mass is increasing and ^ < |, the turning point 
is further and further from the horizon, moving toward the "left" infinity 
(on the other side Einstein- Rosen bridge). And as a limit we get geometry 
of the eternal black hole everywhere except at the left infinity which is now 
singular. 

Another limiting case is a completely closed world with zero total mass 
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Figure 4: 
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Figure 5: 
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Figure 6: 
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Figure 8: 
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Figure 9: 
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m = 0. We have for such a situation 

m = 

(2.14) 

M = 2po/G 

For completeness let us consider the thin shell with negative total mass 
m < (but with positive bare mass M). From Eqn.( pr7D it follows imme- 
diately that aout = — 1 everywhere, the radii thus decrease outside the shell 
up to the singularity at i? = 0. The space-time topology is x rather 
than x R^ as it is for positive total masses m. So, in this case we have 
no spatial infinity at all. In what follows we will consider only positive total 
masses. 

2.3 Black hole universality 

It was already mentioned that the most general black hole is the stationary 
axially symmetric space-time manifold, describing by a very limited set of 
parameters, namely, by its mass, angular momentum and (locally) conserved 
charges. But what does happen when the initial state of the collapsing mat- 
ter is not exactly spherically symmetric (here we consider the case of zero 
angular momentum)? The huge work of many researchers has led to the 
following result (the detailed review of this topic can be found in the book 
|T5|). In the process of black hole formation everything that can be radiated 
is radiating away. It concerns all the multipoles (including the monopole mo- 
ment) of scalar fields, the higher multipoles of electromagnetic, nonabelian 
and gravitational fields. In some theories the scalar monopole can not be ra- 
diated away, but only in those cases when it is rigidly linked to the monopole 
moment (Coulomb part) of the corresponding nonabelian field. This fea- 
ture of process gravitational collapse was first understood by J.A.Wheeler 
and formulated as the (now famous) statement that "The black hole has no 
hair" . Thus, the process of black hole formation is accompanied, in general, 
by the loss of information about the initial state and the preceding history 
of collapsing matter. 

From the very definition of the black hole it is clear that nothing can 
go off the event horizon. But this does not mean automatically that the 
black hole mass cannot be decreased. The matter is that in case of nonzero 
angular momentum there exist around the event horizon a special region, 
called ergosphere, inside which there can be states with negative energy (from 
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the point of view of a distant observer). Then, the following process can take 
place. The particles that fall into the ergosphere, may decay into pairs, and 
one of the component whose energy is negative, may fall into the black hole, 
while another one would go out of the ergosphere (but not out of the black 
hole!) and escape to the infinity, thus decreasing the black hole mass. In the 
literature it is known as the "Penrose process". If instead of particles, we 
are dealing with the fields, then the analogous process also will take place, 
it is called "superradiance". Both Penrose process and superradiance result 
in decreasing the black hole angular momentum. And the real universality 
of black holes is that fact that it is impossible by any classical process to 
decrease the area of the black ole boundary, i.e. for the event horizon area 
A we have always 

dA>0 (2.15) 
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3 Quantum Era of Black Hole Physics 



3.1 Black hole thermodynamics 

In this section (and only in this) subsection we will use some formulas for 
a general axisymmetric black hole with total mass m, electric charge Q and 
angular momentum J. For such a black hole Being the final state of the grav- 
itational collapse, black holes have very interesting properties.lt appeared 
that black holes, during their formation, emit everything that can be emitted. 
Thus, the black hole resembles the body in thermodynamical equilibrium. In 
1972 Jacob Bekenstein suggested that this is not merely a coincidence, and 
the black hole really brings some amount of entropy and has some definite 
temperature. This point of view was supported by the following inequality al- 
ready known at the time, namely, that the area of the black hole horizon can 
not decrease. This resembles the second law of thermodynamics. J. Bekenstein 
suggested that the black hole entropy is proportional to the black hole area 
A (more general, the entropy may be equal to some increasing function of 
the area). 

In this section (and only in this) subsection we will use some formulas for 
a general axisymmetric black hole with total mass m, electric charge Q and 
angular momentum J. For such a black hole 



Moreover, it appeared that the first law of thermodynamics is also valid and 
can be expressed by the following mass formula. 



Here Sm is the mass difference between two stationary black holes with 
slightly different areas SA, angular momenta SJ and electric charges SQ, 
while Q is the black hole's angular velocity and $ is the value of the Coulomb 
potential at its surface. The prefactor k. in front of 5 A is the so called surface 
gravity which equals to 




(3.1) 




(3.2) 



5m = —SA + nSJ + 



(3.3) 



47r 




(3.4) 



1^ = —r 
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Following J.Bekenstein, the surface gravity should be proportional to the 
black hole temperature O, k = aQ. 

Let us consider a Schwarzschild black hole as a thermodynamical sys- 
tem and put it in a thermal bath at some fixed temperature. Then, if the 
temperature of the black hole is the same as that of the thermal bath, the 
amounts of the absorbing and emitting radiation are equal, the black hole 
and the bath are in thermal equilibrium. But this equilibrium is unstable. 
The reason for this is the negative heat capacity of the Schwarzschild black 
hole. Indeed, from the definition of the heat capacity we have 

C = e— = -SnGm^ < (3.5) 

Let us suppose that the black hole is initially in thermal equilibrium with 
the environment. Then, if by some fluctuation the mass of the black hole 
will become a little bit larger than its equilibrium value the temperature 
will become lower than tat of the thermal bath causing the absorption of 
radiation, further increase of the mass and decrease of the temperature, thus 
resulting in the unbound growth of the black hole size and mass. If the initial 
fluctuation results in the small decrease of the mass and, consequently, in the 
increase of the black hole temperature, the black hole starts to emit radiation 
more than to absorb it, thus further decreasing the mass up to the very end. 
In this case the black hole disappears completely. 

In the case of the Reissner- Nordstrom black hole the situation is different. 
The temperature is zero for the extreme black hole (m = \Q\/G), reaches its 
maximum and then decreases to zero when the mass infinitely increases. The 
heat capacity is positive near the extreme state, first grows to infinity, then 
becomes negative hke in the Schwarzschild case. We see that the charged 
nearly extreme black hole is thermally stable like the ordinary bodies. 

3.2 Hawking radiation 

The quantum era of the black hole physics began in 1974 with Hawking's dis- 
covery of the black hole evaporation. Stephen Hawking considered quantum 
vacuum fluctuations of matter flelds on the Schwarzschild black hole back- 
ground and showed that the very existence of the event horizon leads to the 
radiation flow from the black hole (violating the classical theorem dA > 0). 
The surprising fact is that this radiation is nothing more but the blackbody 
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radiation with the temperature 

e^^ (3^6) 

Thus, the black hole entropy is exactly one fourth of the dimensionless hori- 
zon area ^ ^ 

SbH = 7 72" (3-7) 

which is in full agreement with the Beckenstein's conjecture (we use the units 
in which the Planckian constant /I, the light velocity c and the Boltzmann 
constant k are put equal to one; Planckian length equal to Ipi = ^^/hG/c^ = 
G^/^, and Planckian mass is nipi = ^hcjG = 

Remarkably enough the same value of the black hole temperature can 
be obtained in a quite different way. By making a Wick rotation of the 
M2-part of the Schwarzschild or Reissner-Nordstrom line element written in 
the curvature coordinates we obtain a two-dimensional Euclidean half-plane 
(— oo < tE < oo, R+ < R < oo). If we place the origin at i? = assume 
periodicity of the Euclidean time ts and identify points over the period, we 
then get a two-dimensional surface with, in general, a conical singularity 
at the origin. This singularity can be removed by a suitable choice of the 
time period. And the inverse of such a period is equal to the temperature 
in finite temperature quantum field theories and in our case it can be called 
a topological temperature. Practically, it can be done as follows. The line 
element of the two-dimensional Euclidean surface is 



df = Fdr' + Idr^ 
F 



2Gm Ge^ I , 

1 - —;r- + 7^ = ;:2 - ^+)(^ - ^-)' 



(3.8) 



(3.9) 



e2 



r± = G \m±\lm? - —\ . (3.10) 

The Euclidean time r is assumed to be a cyclic coordinate. In the vicinity 
of the horizon F = (r — r+)(r — r_)/r^ . Introducing a new coordinate 
p — 2r+y'r — r+/-^r+ — r_, we get 

df = ^l±^lzjLp^dT^ + dp" = dp" + p^dif^ (3.11) 
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And this is nothing more but the metric of a locally flat two-dimensional 
surface, written in the cylindrical coordinates. The point p = is, in general, 
a conical singularity. To remove it, we have to require the period of the 
azimuthal angle (p to be 27r. Thus, 

r = 3.12 

r+ — r_ 

Its inverse, 

2G'./m2 - 

e = — L ^ (3.13) 

equals exactly the temperature of the Reissner- Nordstrom black hole. 

All this means that quantum matter fields in the black hole background 
arc described, in fact, by the finite temperature quantum field theory. In this 
sense, the appearance of the black hole temperature looks quite naturally. 

The violating of the classical law > does not mean that the radi- 
ation can be emitted off the black hole, because the black hole itself is still 
considered as the classical object. The calculations of the stress-energy ten- 
sor of the vacuum fluctuations show that their energy density is negative. 
Therefore, the negative energy flow falls into black hole, thus diminishing its 
mass, while the, equal to it, positive energy flow escapes to infinity. One can 
consider such a situation as that the quantum fiuctuations of matter fields 
form something like an ergosphere in the vicinity of the event horizon, and 
inside this ergosphere the Penrose process (for particles) or the superradiance 
effect (for radiating fields) may take place. 



3.3 Rindler's space-time. 

Except the black holes there is one more example of the space-time manifold 
which observers "see" a black body radiation. This is the Rindler's space- 
time. It is locally ffat, therefore, sufficiently simple and all the quantum 
field theory equations can be exphcitly solved and investigated in details. 
On the other hand, this example is very important since there is exists a 
deep relation between Rindler's and black hole manifolds. Let us consider 
the two-dimensional Minkowski space-time with the metric 
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The Rindler's coordinates are obtained by the following transformation 

t = ±—6""^ sinh ar] 
a 

(3.15) 

a; = ±-e"^ cosha?7, (3.16) 
a 

The upper (lower) sign is for x > (<)0. The line element now reads as follows 

ds^ = e^^^dr]^ - d^^) (3.17) 

Rindler's observer, sitting ai ^ = const., experience a constant proper accel- 
eration a = ae~"^. Being locally fiat, the Rindler's manifold differs, neverthe- 
less, from the Minkowski manifold, it has different boundaries. The boundary 
of the Minkowski space time are trivial infinities (spatial, temporal and null). 
Rindler's coordinates cover only one half of the Minkowski space-time, cor- 
responding to |a;| > \t\, and the lines x = ±t are the event horizons that 
serve as new boundaries (in addition to the infinities). Exact quantum field 
theory calculations show (the details can be found in the book |]TB[) that the 



Rindler's vacuum is, actually, with the black body radiation which tempera- 
ture equals the so called Unruh's temperature Tq = Every local observer 
measures its own local temperature T = To(5'oo)~^^^ = Tqc'"-^ = 

It should be stressed that the appearance of the temperature (and heat 
bath) in Rindler's space-time is a global effect, depending essentially on the 
boundary condition at the event horizon. And we can do here the same trick 
as we did with the black hole space-time. Namely, we can introduce the 
cyclic imaginary time rj irj = t and spatial coordinate ( = -e""^ the line 
Rindler's element now becomes 

dl^ = a\V + dC. (3.18) 

The obtained two-dimensional euclidean surface has, in general, a conical 
singularity that disappeared only if the period of our imaginary time r equals 

The inverse to this period is exactly the Unruh's temperature. Note,that 
procedure, described above, is of global character, because the choice of the 
period in the point dictates actually the period of the "proper" imaginary 
time for all Rindler's observers. 

Let us consider now a static spherically symmetric space-time the line 
element of which can always be written in the form 

ds^ = e^dt^ - e^dr^ - r^dQ^, (3. 19) 
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where v and A are function of the radius r. By the equivalence principle, the 
static observer experiences a constant acceleration. And indeed, it can be 
easily shown by direct calculation that the acceleration equals a = Ke~'^^'^. 
where 



^e— (3.20) 

is the so called surface gravity. A static observer, sitting at r = ro 7^ can 
forget about the angular part of metric ( |3.19|) and consider a 2- dimensional 
static curved space-time, such a 2-dimensional surface can be easily trans- 
formed to the local Rindler-like coordinates. Indeed, 

ds^ = e^dt^ - e^dr^ = e%dt^ - e^-^dr^) = e^'idrf - d^^), (3.21) 

where t] = t and d^ = . We can calculate the acceleration parameter a 
in the following way. Let r = rg + 6r, then 6^ = e^~6r The conformal factor 
near r = tq equals e^(ro)(l + u'Sr) = e^"^°(l + 2aS^). From this it follows 
that e'^°'^° = c'^Itq) and a = ■Ye^~(ro) which is exactly the surface gravity k 

The equivalence of the Rindler's observers to the static ones in the spher- 
ically symmetric space-times does not mean that the latter is heated. This 
equivalence is only local. To have the global equivalence (of course, with 
different Rindler's observers for different static observers) we need the event 
horizon, and that is exactly what we have in the case of spherically sym- 
metric black holes. But this is not the end of the story. The reduction 
of the four-dimensional spherically symmetric space-time to 2-dimensional 
{t, r)-surface can be safely done only if r 7^ 0. The point r = is the coor- 
dinate singularity. Thus, the equivalence described above is valid only if the 
corresponding euclidean surface does not contain such a point. Otherwise 
we would have to consider spherically symmetric Rindler's manifold, which 
essentially non-static. 

And, finally, the last note. The temperature (if it exist) measured by a 
local static observer obeys the law Q^/goo = const. Therefore, the temper- 
ature Go = ^ is measured by observer, sitting at the point where goo = 1- 
For the asymptotically flat space-time this is, of course, the point at spatial 
infinity. 
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3.4 Quantum mass spectrum 

Historically, the discovery of the quantum nature of radiation has led to 
the construction of the hydrogen atom model and to the explanation of its 
discrete energy spectrum. Analogously, the quantum nature of radiation 
suggests that the mass spectrum of black holes is, actually, quantized. Jacob 
Bekenstein proposed the equidistant spectrum for the black hole horizon 
area. His motivation is as follows. The entropy value is just the logarithm 
of the number of ways in which we can construct the black hole with given 
parameters and related to the amount of information hidden inside the black 
hole. The amount of information is naturally quantized. Therefore, the black 
hole area which is proportional to the entropy should have an equidistant 
spectrum. A more quantitative explanation is due to Slava Mukhanov. We 
generalized his result originally derived for the Schwarzschild black hole, to 
the Reissner-Nordstrom charged black hole. Let the black hole emits the 
quantum energy uj due to transition from the n-th state to the {n — l)-th 
state. The "typical" energy of this quantum minus the work done against 
the Coulomb attraction should be proportional to the temperature, i.e. 

0Jn,n-i - <^dQ = Trae, (3.22) 

where dQ is the charge of the quantum in question, and a is a slowly varying 
function of the integer quantum number n. Inserting, then, Eqn.( p^22D into 
Eqn. (|3.3|) we have {(jJn,n-i = dm) 

dA = AirGa, (3.23) 

or 

dA = AnGadn, (3.24) 

Integrating we get 

A = ATiGan + AnGC, (3.25) 

where C is a constant of integration, and a is some another slowly varying 
function of n. It is natural to suppose that the value n = corresponds to 
the minimal possible value of the horizon area, e.i., it should correspond to 
the extreme black hole. This gives us G = Q^. It is clear that the Hawking's 
formulae are valid only if the back reaction of the radiation on the space- 
time metric can be neglected, e.i., for law (comparing to the black hole mass) 
temperature. Thus, we are interested in two different cases, in the large black 
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hole regime and in the nearly extreme regime. In the first case and 
we assume the following expansion 

a = ai + /3i^+7i^ + ..., (3.26) 

— «i, 

n 



while in the second case n <^ Q"^ and we have 

n n' 



a = a2 + /327^ + 727^ + ---, (3.27) 



— > 1. 

n 

The corresponding expansions for the black hole mass are 
m = — -i=\joi.\ I 1 + 



2a/G V 2q;i n 

71 (/3i + i)(/ii + 4) \ 

2q;i Saf / 



+ ^ ' \^^---\^ (3-28) 



2 



and 



{ 2 + — 7^ + --- ) , (3.29) 



2Vg V 4 



n< g^ 



The values of a's, /3's and 7's can and will be determined in the next Section 
for the specific model. 



4 Quantum Black Hole Models 
4.1 Eternal black hole 

We have already described in detail the classical eternal black hole and men- 
tioned that it has no dynamical degrees of freedom. In this respect it is anal- 
ogous to the Coulomb field in electrodynamics. Nevertheless it is instructive 
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to quantize such a system. Of course, the full program of quantization can 
not be fulfilled because quantum gravity is not renormalizable pertubatively. 
But the quantization of the vacuum spherically symmetric space-times can 
be done in a nonperturbative way in the so-called frozen formalism. In this 
formalism we quantize only one, radial, degree of freedom. Though not com- 
plete, such a procedure can give us rather reliable results, especially in the 



light of the recent result |T8| which support the suggestion that the four- 
dimensional Einstein gravity may appear nonperturbatively renormalisable 
|19| . Quantization of the eternal black hole in frozen formalism was done by 



H.Kastrup and T.Thiemann |^ in the so-called quantum loop gravity and 
by K.Kuchar who used the conventional geometrodynamics. Even for 
such a simple space-time manifold the problem appeared nontrivial. But the 
final result was the same in both approaches, namely, the quantum functional 
on spherically symmetric vacuum metrics depends only on the Schwarzschild 
mass. And in the absence of dynamical degrees of freedom it is impossible 
to extract any information about a black hole spectrum. Physically such a 
trivial result can be understood if one compares the eternal black hole with 
its electro-dynamical analog, the hydrogen atom. Let us imagine that we first 
allow the electron to "fall" classically on the proton and then start quanti- 
zation procedure. What we would get? Absolutely nothing. Of course, in 
electrodynamics we have both positive and negative charges, and the mass 
of the eternal black hole can only be positive. And, instead of "nothing" the 
curvature singularity is developing in the process of gravitational collapse. 
But the qualitatively, the picture is the same. Therefore in order to obtain 
physically meaningful results we have to "switch on" some dynamical de- 
grees of freedom. In subsequent Subsection we consider the quantization of 
spherically symmetric thin dust shells, for which we already developed the 
classical theory. 

4.2 Naive quantization 



The material of this section is based on the works l^^j, ^ ^|. Now, we 



investigate a specific quantum black hole model, which classical counterpart 
is described in the first part of the lectures. It is a self-gravitating spheri- 
cally symmetric thin dust shell with bare mass M, total mass m and electric 
charge e. The "correct" quantization procedure (which starts from the "first 
principles" ) we postpone till the next Chapter and consider here a very sim- 
ple quantum model based on the classical equation of motion for our thin 
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shell. Since m is a total mass of the system and the Eqn. ( |2.8| ) represents the 
energy conservation law it is reasonable to consider it equal numerically to 
the Hamiltonian of the system in question. Thus, we have 

m = M^p^ + 1 (4.1) 

2p 

To quantize a system we should construct first the classical Hamiltonian 
and then to replace the classical Poisson brackets by the quantum commu- 
tation relations. Given the classical expression E{p,p) for the energy of the 
system as a function of its radius and velocity (actually, rapidity in our case), 
the conjugate momentum p, Lagrangian L and Hamiltonian H can be found 
from the following relations 

dL 

E = pp-L = P—-L, (4.3) 

which give for the conjugate momentum 

f dEdp 

Substituting for E the Eqn.([4.lD one gets 

p = M log(p + + 1) + F{p), (4.6) 



L = M(plog(p + + 1) - Vp' + 1) + pnp), (4.7) 

where F{p) is an arbitrary function. The choice of this function does not af- 
fect the Lagrangian equations of motion, and leads to the canonically equiv- 
alent systems in the Hamiltonian approach. This choice may become impor- 
tant when one explores the rather complex structure of the maximal analyt- 
ical extensions of black hole space-times trying to construct a wave function 
covering the whole manifold. But here we do not need such a complication 
and put F{p) = 0. From Eqn.([4.6[) we get for p 

p = sinh^, (4.8) 
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and, after inserting this into Eqn.( |4.1| ) we obtain the following Hamiltonian 
H. 

H = Mco^\i47-^ -■ (4.9) 

It is convenient to make a canonical transformation to the new variables (M 
is the bare mass of the shell) 

x = Mp, n = ^p, (4.10) 

then 

ij' = M(coshn- 2x 1' 
Let us now impose the quantum commutation relation 

[U,x] = -i. (4.12) 
Then, using the coordinate representation with 11 = —id/dx and the relation 

e-^J^^(a;) = ^(x -i) (4.13) 
we obtain the Schroedinger stationary equation 

H^{x) = E^{x), (4.14) 



M ( GW^ - \ 

— {^{x + z) + ^{x -i) ^(x) = rn^ix). (4.15) 

2 V X ) 

Introducing the notation e = m/ M we can rewrite this as 

^{x + i) + ^!{x -i) = \2t + j ^(x). (4.16) 

Some remarks are in order. First, the Schroedinger equation obtained is 
not the differential, but the equation in finite differences. This is a direct 
consequence of the quantization in proper time which is quite natural in the 
framework of General Relativity. Second, the total mass m is not arbitrary 
anymore but it is an eigenvalue of the Hamiltonian operator, subject to the 
condition that for corresponding eigenfunctions the Hamiltonian is Hermitian 
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on the positive semi-axes. Third, it should be noted also that the step in our 
finite difference equation is along the imaginary axes, so the solutions should 
have certain analyticity to be continued to the complex plane. 

The Hamiltonian picture we got is not equivalent to that which could be 
obtained directly from the Einstein-Hilbert action. But our procedure self- 
consistent. We started from equations of motion written already in proper 
time and constructed the new Lagrangian which results in just the equations 
in proper time with no restrictions or/and constraints. We do not expect that 
the resulting mass spectrum will be correct. But our model is extremely sim- 
ple and, as we shall see soon, is exactly solvable, and it seems very instructive 
to develop here some new methods and ideas. 

We consider the non-relativistic limit of the obtained Schroedinger equa- 
tion. With all dimensional quantities restored the shifted argument [x ± i) 
becomes 

x±i-^ M ( p± -J-z ) ^ p± ^ ■ 



M J ' Mc 

so for p ^ h/Mc we obtain the non-relativistic limit. Expanding Eqn.( ^?T6D 
up to the second order in {h/Mc) we get 

^ = (E-M)^ (4.17) 

which is just the non-relativistic Schroedinger equation for the radial s-wave 
function, {E — M) being the non-relativistic energy of the system. For the 
negative values of {E — M) acceptable solutions exist only for a discrete 
spectrum of energies, namely 

M(GM^ - 6^)2 

{E-M)^ = -^^^^^^^^^^, n=l,2,... (4.18) 

For G = it reduces to the well known Rydberg formula for the hydrogen 
atom. 

Now, we show how to solve the obtained Schroedinger equation in finite 
differences. And we will do this first in the momentum representation. 
It is convenient to introduce the following parameters, A and (3, 

e = cosA, a = GM^ - = 2/3sinA, 

in terms of which our Schroedinger equation becomes 

^(x + i) + ^{x-i) = 2 (cos A + ] ^{x). (4.19) 
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Since for the black holes a > 0, the signs of A and beta are the same. We 
choose A > and, hence, /? > 0. Let us return to the operator form of the 
equation, 

(cosh p - /5 sin Ai:"-^) \& = cos A\&. (4.20) 
In the momentum representation operators p and x act as follows 

d 

p^p = p^'p, x^'p = i-Q^'^p, (4.21) 

where I Psip is a wave function in the momentum representation. The operator 
x"^ is not well defined. The ambiguity can be removed by adding suitable 
terms to the potential which are proportional to (5-function and its derivatives 
at the origin. But, instead, we can multiply the equation by operator x from 
the left. By doing this we get 

d 

i— (coshp - cos A)^p = /3sin A^p, (4.22) 

or 

d sinhp + i/3sinA 

— log ^p = — . 4.23 

op cosh p — cos A 

After introducing a new variable, z = e^, the Eqn.(^]2^) takes the form 

9 , , z^ + 2iPsmXzi 1 p+1 p~l 
-TT-log^ip = — — — = + -, 

OZ Z[Z'^ — 1 COS Xz ^\) Z Z — Zq Z — Zq 

where zq = e*^, Zq = e"*^. (4.24) 
The above equation can be easily solved, the result is 

^ \ /3 



{z - Zo){z - Zo) \Z - ZoJ 

z = eP. 

This solution has very important property, it is periodical with the pure 
imaginary period 27ri. This property will be explored below. 

Now we will transform the solution found in the momentum representa- 
tion, to the coordinate representation. It can be done by the inverse Fourier 
transform, 



^{x) = —= I e'P''^pdp. (4.26) 
v27r 
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We obtained before the unique (up to the multiphcative constant) solution 
for , namely, Eqn.(^]26|). But due to the periodicity of \l/p we can shift 
the argument p — > p + 27rki{k = 0, ±1, ±2, ...) which will result in the shift of 
the path of integration in the complex momentum plane from the real axis 
to the parallel one. And after such a shift we again obtain a solution in the 
coordinate representation. But, 

oo 

^i(x) = [ e'P'^'^pdp, 

\/2n J 



oo 

vl>2(x) = -^ / e'^P+'^'''^^^p+2.Hdp = 
V 27r J 

— oo 

oo 



-oo 



^2(x) = e-'"'="^i(x). (4.27) 

Thus, given one solution (say, \l/o), "we can construct in this way a countable 
number of solutions solutions, and, in general, 

"^.eueralix) = f Cfce'^-'^- j <ilo{x). (4.28) 

\fc=— oo / 

The infinite sum in parenthesis is nothing more but a Fourier series for a 
periodical function with an imaginary unit period i. And this reproduces the 
following property of the solutions in coordinate representation. Given some 
solution, say, \E'o we can construct new solution, multiplying it by any 
periodical function C{x) with the imaginary period i. 

How many different solutions do we need to construct the general solu- 
tion? The very fact that in the momentum representation we obtained es- 
sentially one solution proves that for this we need only one solution \I^o which 
we call the fundamental solution (though it deserves the name superfunda- 
mental). Now we will construct the particular fundamental solution by the 
suitable choice of the contour of integration in the inverse Fourier transform 
integral. Note, first of all, that our solution, Eqn.( ^l26D , in momentum repre- 



sentation has countable number on branching points in the complex momen- 
tum plane which can be combined in pairs, {iX + 2TTki, — zA + 27r(/c + k = 
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0, ±1, ±2, .... Connecting two branching points in each pair by a cut we ob- 
tain the complex plane with countable number of cuts. On the corresponding 
Riemann surface our solution is a single valued analytical function of com- 
plex variable. Our choice of the contour of integration is as follows. We 
will integrate first along the real axes from left to right (i.e., from — oo to 
+oo), then along a short curve at the right infinity {p p + 27ii,p —>■ +oo), 
then along the straight line y = 2TTi, parallel to the real axis, from left to 
right, and finally along a short curve at the left infinity back to the real axis. 
The integration along the short curves at infinities between the straight lines 
gives zero contributions for positive values of x (for negative x we can choose 
the straight line y = —2Tii instead oi y = 27ri). The integration along each 
straight line gives us a solution the linear combination of which is again a 
solution. Thus, the inverse Fourier integral along such a closed contour gives 
us a solution. This contour can be distorted to become a contour around the 
cut (Ai, {2ti — X)i) for a; > (or around the cut — (27r — for a; < 0). 

Thus, 

■^o{x) = (f e'f"''^pdp, x>0, (4.29) 



where 



2n 

C- 

= (p e'P^'^Updp, X < 0, (4.30) 



^2n 

c. 



% = C- t(^)' (4.31) 

[z - zo){z - Zo) \Z-ZqJ 
z = eP, Zo = e'\ Zo = e"*^ 

In what follows we restrict ourselves to the case of positive x only. The above 
integral representation for \E'o can be simplified if we integrate Eqn. (|4.29|) by 
parts. Using the fact that 

' dL,^), (4.32) 



{z — zo){z — Zo) 2isinA \ z — zo 
we get 

c+ 
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(the extra term vanishes because our contour has no boundary and for conve- 
nience we omitted some constant factors). We will use use Eqn. ([4.33|) as an 
integral representation for the fundamental solution of our finite differences 
equation, Eqn. (|4.19|) . 



Let us investigate the fundamental solution in more details. Our aim 
is to reduce the integration along the closed contour C+ around the cut 
to the finite interval between the corresponding branching points. But for 
/3 > 1 (remember that we have chosen (3 positive) we have the non-integrable 
singularity at the lower integration limit. To avoid this difficulty we need 
some recurrent relation for lowering the parameter /3. To find such a relation 
we integrate by parts Eqn.( 4733D , the result is the following. 



mp{x) = ^p-i{x) + j-^ {^6eta-i(x) - e-'^^^_i(x + l)] . (4.34) 

From the above relation it is easy to derive the structure of for general 
values of /? > 1. Let us take /? = n + /3, with /? < 1. Then, ^ p \s the sum of 
two terms, which of them is the product of some polynomial of n-th degree 
and of ^{x) or \l/^(x + i). Thus, we can proceed assuming 13 < \. First of 
all, consider separately the case (3 = 1. We have 

•qf^ = X (p 6*^""— —dp = 2mx hm ^ j 



e'P - e*^ p-^ix eP - e'^ 



a 



+ 



= - {Ane-'^ sin X) xe-^"" . (4.35) 
From this it follows that all for positive integer (3 = n has the form 

^„ = P„(x)e-"^ (4.36) 
where Pn{x) are some polynomials of n-th degree. In the case < /? < 1 



the integrand in the Eqn.( 4.33 ) is integrable on both ends of the cut in the 
complex momentum plane and we are able to convert the integral along the 
closed contour C+ into the integral along the finite interval. We get the 
following, 

^p{x) = j e'^^Udp = (1 - e^-^) x^p{x), (4.37) 



where 



i{2-K-X) 

eP - e 



$^(x) = / e*P^ ( — ^ ] dp, x> 0. (4.38) 



iX 
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Changing the variables, 

e'i -e'^ = -{2ismX)y, (4.39) 

we get 

1 p 

= (1 - e-'^') e-''^ /(!-(!- ^'"'') vY^'' (^) dy. (4.40) 



Comparing this with the well known integral representation for the Gauss's 
hypergeometric function, F{a, b; c; z), 

1 



FMc;z) = t'~\l-tr-'-'{l-tz)-'^dt, 



1 

B(.,.)^/*-(l-*r'*, (4.41) 



we see that the integral in Eqn.(OO) is the Gauss's hypergeometric function 



F(a, b; c; z) with the following values of parameters, 

a = l-ix, 6 = 1-/3; c = 2; z = l- e'^'^, (4.42) 

and 

5(6, c - 6) = 5(1 - /5, 1 + /5) = (4.43) 

sm vrp 

Finally, for the fundamental solution i3{x) we get 

^^(a;) = (-47r/5e-^^ sin A) a;e-^^F (l - zx, 1 - /3; 2; 1 - e'^^^) . (4.44) 

The above expression was derived for < /5, 1 only, but the same is valid for 
any value /5 > 0. And, finally, it can be easily shown that the original equa- 
tion in finite differences, Eqn.( ^rT9D , is a direct consequence of the following 



Gauss's recurrent relation 

{2a—c—az+bz)F{a, 6; c; z) + {c—a)F{a—l, 6; c; z)+a{z—l)F{a+l, 6; c]z) =0 

So, we found the general solution to the equation in finite differences with 
Coulomb potential, 

^(x + + ^(x - z) = 2 l^cos A + - sin A^ ^(x). (4.45) 
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Now we should remember that this equation is actually the radial Schroedinger 
equation for zero angular momentum, and the coordinate x runs from zero 
to infinity, i.e., x takes values in the positive semi-axis only. Thus the Hamil- 
tonian, 

should be the selfadjoint operator on the positive semi-axis rather than on 
the whole real axis as dictated by the quantum mechanics postulates. In 
addition, the wave function should be square integrable on the semi-axis. 
The corresponding extension of the above Hamiltonian was found by Petr 
Hajicek. It appeared the wave function should obey the following boundary 
conditions at the origin, x = , 

^(2«)(o) = 0, n = 0,l,.... (4.47) 

That is, the function itself and all its even derivatives should be zero at the 
origin. The appearance of infinite number of conditions is due to the infinite 
order of the operator. P. Hajicek found also a conserved probability current 
J{x) for the Schroedinger equation in finite differences, Eqn.( f4.45| ) . It is 
given as usual by the equation 

J'(x) = i (^*H'^ - ^^H^' 

so that 

(^*^)- + J' = Q 



Writing the Hamiltonian ([4.46[) as follows 



H = lytll^m^^^_^ (4.48) 
2^ (2k)\ ^ ^ 2x ^ ^ 



we obtain for J(x) 



oo / 1 \^ k 

j(^) = , [^*^'~'^^(2A; - /) - . (4.49) 
fc=0 ^ i=\ 

P. Hajicek showed that the boundary condition ( ^.47| ) implies that 

J(0) = 0. (4.50) 
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The boundary condition ( [4.47 ) is a direct generalization of the boundary con- 



dition in the case of non-relativistic Schroedinger equation which is of second 
order. In this particular case the boundary conditions ( |4.47] ) are reduced to 



the single one, \E'(0) = , which is together with the square integrability 
condition ensures existence of a discrete spectrum for bound states. And, in 
general, if we expand Hamiltonian ( 4.46| ) in series of derivatives and consider 



the truncated Hamiltonian of 2A^th order, the corresponding Schroedinger 
equation becomes the differential equation of 2iVth order. It can be shown 
that the asymptotics at of such an equation at the origin are 

^ ~x,a;^...,a;2^-^l + x2^-Mogx (4.51) 

and at the infinity they are of the form 

where are real for bound states. The general solution is defined up to a 
normalization factor and depends, actually, on {2N — 1) arbitrary constants 
which are to be determined using boundary conditions. The square integra- 
bility condition reduces the number of unknown constants to (A^ — 1). But 
for the truncated operator of 2iVth order we have N conditions at the origin. 
They can only be satisfied for a discrete spectrum of eigenvalues. Note that 
the last of the asymptotics ([4. 511) is excluded by the boundary condition 



at the origin. We can reverse the above consideration and start from the 
asymptotic behavior at the origin. To satisfy the boundary conditions there 
we have to exclude the last asymptotics ( |4.51| ) and all the even terms in the 



expansion up to the (2A^ — 2)th order. By doing this can determine N of 
(2A^ — 1) unknown constants. Thus we are left with only {N — 1) constants 
for N boundary conditions at infinity. But the situation is not so simple in 
the case of an infinite order operator. First of all, how to get rid of the loga- 
rithmic term in asymptotics? The problem is that if we first go to the infinite 
order limit, then the logarithmic term unconditionally disappears (|x| < 1!) 
leaving 1 as an asymptotics. But the latter can be easily reproduced by some 
infinite linear superposition of the remained asymptotics. Subtracting one 
from another we would obtain the solution which would satisfy formally all 
the boundary conditions at the origin still having enough arbitrary constants 
to satisfy boundary conditions at infinity. But if we first differentiate 2A^ 
times the asymptotic solution containing logarithmic term we will get (l/x) 
term. Thus, requiring that the "good" solution should be infinitely differen- 
tiable we can reach our goal. Moreover, we have to require the analyticity 
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of the wave function at least on the real axis. Thus is because we implicitly 
have used the analyticity of the solutions in transition from the differential 
equation of infinite order to the finite differences equation. The situation 
is not good at the infinity either. We saw that in the case of our partic- 
ular equation we can obtain a new solution multiplying the known one by 
exp {—2'Kkx) . Choosing the large enough value of k we are able to convert 
the "bad" at infinity solution to the "good" one. We will see in a moment 
that it is the analyticity condition which solves this problem. 

Let us consider the asymptotics of the general solution at x — and 
X — s> oo . The idea is the following. First, we find the asymptotical behavior 
of the fundamental solution at x — > and x — oo. We have, 

x\ r = l,2.... (4.53) 

at X — > and 

X ^ oo (4.54) 

It is interesting to note that the same asymptotics can be found directly from 
the original Schroedinger equation in finite differences without knowing the 
exact solution. 

Comparing the asymptotics at x = and at infinity we see that the 
fundamental solution is an analytic function only if 

(3 = n (4.55) 

(remember that we chose /5 > 0). This leads to the discrete spectrum for the 
total mass m. Indeed,we have 

13 = = n, a = kM^, 

2 sin A 

m 

e = — = cos A, 
M 



1 - 



For (3 = n the hypergeometric series in Eqn ([4.44|) terminates and we are left 
with the following fundamental solution, 

^„(x) = P„(x)e-"^ (4.57) 
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where P„'s are some polynomials of order n. They can be calculated directly 
from the definition, or from the generating function 



n=0 ■ n=0 



V-2 + cotA-«y ^ ^ 

At the end of this section we give an example of solving the boundary 
conditions and find the wave function of the ground state (i.e., for n = 1). 
The fundamental solution for n = 1 was found to be 

= xe"^^ (4.59) 

(we omitted here the irrelevant constant factor). The general solution obey- 
ing the boundary condition at infinity (exponential falloff at x — » oo) has the 
form 

oo 

^igen = xe-^^ J]cfce-2-^^ (4.60) 

and the coefficients are to be determined by solving boundary conditions 
at the origin, 

^(20(0) = 0, / = 0,1,... (4.61) 

Differentiating Eqn.( |4.60|) 21 times we get the following infinite set of linear 
equations, 

oo 

2l^Ck{X + 2TTky^-^ = 0, 1 = 0,1,... (4.62) 

A:=0 

The determinant of this system is identically zero because the first line con- 
sists of zeros. We can truncate this system at some specific value of /, calcu- 
late all the determinants and minors and then, after appropriate renormal- 
ization (say, putting cq = 1), take the limit / — > oo. In our simplest case of 
the ground state it is possible to make all the calculations up to the very end 
with the following result 

(_l)fer(^ + A;) 
fc/ ^ = 0,1,... (4.63) 
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Inserting this into Eqn. ( ^^BOD we can write the ground state wave function 
(up to the normahzation factor) in a very simple and nice form, 



fc=0 



k=o \ / \ 



xe 



2~^--^-T. (4.64) 

(cosh TTX) T 

It is easy to see that after shifting the parameter A — >• A + 27r/, where / 
is a positive integer, we again obtain a solution satisfying all the boundary 
conditions (technically such a solution can be obtained if put zero first I 
coefficients Cm, m = 0, ...(/—I) in the infinite set of linear equations considered 
above and normalize to unity q). Moreover, the introduced earlier conserved 
current is identically zero for any superposition (with complex coefficients) of 
these functions. This means that each eigenvalue of a total mass (energy) is 
infinitely degenerate and the appearance of the infinite number of branching 
point refiects this fact. 



4.3 Quantum black holes and Hawking's radiation. 

In the preceding sections we have got a complete solution to the Schroedinger 
equation in finite differences describing a quantum mechanical behavior of the 
self-gravitating electrically charged spherically symmetric dust shell. For the 
bound states we found the discrete spectrum of mass which coincides with 
Dirac spectrum if we put zero the so called the so called radial quantum 
number. Similarly to the case of hydrogen atom, to which our spectrum 
reduces in the non-relativistic limit, we may conclude that the shell does 
not collapse without radiation. The starting classical situation is however 
different. The classical hydrogen atom collapses radiating continuously while 
the classical spherically symmetric shell collapses without radiation. But 
without such a collapse the event horizon and thus the black hole can not be 
formed at all. Moreover, if we calculate the mean value of radius of the shell 
using the wave function we will find that at least for large values of principal 
quantum number it is far away of classical event horizon, so the shell spends 
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most of its "lifetime" outside the classical black hole. So, to get a black hole 
solution we need some new criterion and more detailed investigation of the 
spectrum is needed. Let us consider the dependence of a total mass m on a 
bare mass M , the other variables, e and n fixed. We see that there are two 
branches, an increasing one for 

(GM' - e') (SGM^ - e^) < An^ (4.65) 

and a decreasing one in the opposite case. Note that now 



m 2 

and this value is greater than the corresponding classical value. The increas- 
ing branch corresponds to the "black hole case" while the decreasing branch 
- to the "wormhole case". Using a "quasiclassical" argumentation we can 
say that for the states obeying inequality ( 4.65| ) the expectation value of the 



radius lies outside the event horizon on the "our" side of the Einstein-Rosen 
bridge (thus replacing the notion of "classical turning point" by that of "ra- 
dius expectation value"). For the decreasing branch the same occurs on the 
"other" side of the Einstein- Rosen bridge. It is clear now why the value of 
total mass-bare mass ratio in quantum case is greater than the corresponding 
classical value. The origin of this is just the replacing of the classical turning 
point by the mean value which is smaller thus giving rise to the increase in 
the total mass. Following such a line of reasoning we should conjecture that 
the maximal possible value of a total mass m for charge e and principal quan- 
tum number n fixed corresponds to the situation when the mean value of the 
radius of the shell lies at the event horizon making its collapse possible. The 
further increase in the bare mass m will lead to the wormhole with smaller 
total mass. Generalizing, we can introduce the notion of the quantum black 
hole states. For given values of electrical charge e and principal quantum 
number n the quantum black hole state is the state with maximal possible 
total mass m. Thus, for the quantum black hole states we have, instead of 
inequality ( [4. 651) the following equation 



3{GMY - AGM'^e^ + - An^ = 0, (4.67) 
the solution of which subject to the condition GM"^ — > is 



GM^ = -e^ + -Ve^ + 12^2. (4.68) 
3 3 
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For the black hole mass spectrum we get eventually 



rriBH 



1 IaVS f 2 



6 



n = 1,2,3,... 

For the uncharged black hole with e = we have immediately 

2 

m = ^ \fnmpi. (4.70) 
27 



Two limiting cases are of interest. In the first case 

«1, 



12n2 



and the charge of the black hole gives rise only to the small corrections 
to the mass spectrum of the Schwarzschild black hole. The corresponding 
expressions for the bare mass M and the ratio m/M are, respectively, 

V 3 I 6 n 48 n2 



m rz I a/Sc^ 5 i 



In the second case 



12n2 



»1, 



lei / \ 
The corresponding expressions for M and m/M are, respectively. 



e 
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It is interesting that in this second hmiting case the black hole mass takes 
values nearly the mass of a classical extreme Reissner- Nordstrom black hole. 
It is easily seen that in both limiting cases | A m|/m <^ 1, where | A m| is 
the difference in masses for nearby energy (mass) levels. So, these limits are 
essentially quasiclassical ones, and the corresponding mass spectra should 
be compatible with the existence of Hawking's radiation. This is indeed the 
case. Let us write once more the spectra that follow from the black hole 
thermodynamics in the same limiting cases. 

m = — i=Jai 1 + 



2VG V 2q;i n 
7^ _ (A + l)(A+4) \ 
2ai 8al J 



+ + (4.75) 



and 



lei / an \ , . ^ 



Comparing this with the expansion obtained from the black hole spec- 
trum, we have 

16 ^ 1 5^3 

a2 = 2, /32 = 2, 72 = -15,... (4.77) 

Comparing values for ol\ and 0^2 we see that a is indeed a slowly varying func- 
tion. Thus, we see that in the limiting case of low temperature our black hole 
spectrum is compatible with the existence of Hawking's radiation. Finally, 
we would like to note that our black hole spectrum obeys also the third law 
of thermodynamics. Indeed, for a nearly extreme Reissner-Nordstrom black 
hole the lowest energy level (for n — always exceeds the critical value 
(= |e|/v^), so the zero temperature state is not accessible. 
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5 Geometrodynamics for Black Holes and Worm- 
holes. 



This Chapter is devoted to the "correct" quantization of the spherically sym- 
metric gravity with thin shells. It means that we start with the Einstein- 
Hilbert action reduced to the spherically symmetric space-time, make use of 
the Arnowitt-Deser-Misner (ADM) formalism in order to obtain the canoni- 
cal constraints and Hamilton equations of motion and formulate the quantum 
theory. 

Our model remains the same. This is just a self-gravitating spherically 
symmetric dust thin shell, endowed with a bare mass M. The whole space- 
time is divided into three different regions: the inner part ( I ), the outer 
part ( II ) containing no matter fields separated by thin layer III, containing 
the dust matter of the shell. The general metric of a spherically symmetric 
space-time is now convenient to write in the form: 

ds'' = -N'^dt'' + L\dr + N'dtf + R^de'' + sin^ ^d^^) (5.1) 

where {t,r,9,(j)) are space-time coordinates, N, , L, R are some functions 
of t and r only. Trajectory of the thin shell is some 3-dimensional surface E 
in space-time given by some function r{t): = {(i, r, ^, 0) : r = r{t)}. In 
region I r < f — e, in a region II r > f + e, region III is a thin layer f — e < 
r < f + e. We require that metric coefficients A^, A^'', L and R are continuous 
functions but jump discontinuities could appear in their derivatives at the 
points of E when the limit e — > is taken. The action functional for the 
system of spherically symmetric gravitational field and the thin shell is 

S = Sgr + S shell — / ^'^'^ R\J —gd'^x + (surface terms) — 

-M Jdr ^ ^ 

s 

It consists of the standard Einstein-Hilbert action for the gravitational field 
and matter part of the action describes a thin shell of dust. The surface 
terms in the gravitational action and the falloff behavior of the metric and its 
derivatives were studied in details by Karel Kuchar. So we will not consider 
this question and will use the results of Kuchar when needed. We will be 
interested in the behavior of the action and constraints on the surface E^ 
representing the shell's trajectory. The complete set of degrees of freedom of 
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our system consists of the set of N{r, t), N^{r, t), L{r, t), R{r, t) which describe 
gravitational field and f(t) which describes the motion of the shell. The 
metric ( [5.1| ) has the standard ADM form for 3+1 decomposition of a space- 
time with lapse function N , shift vector = (iV^, 0, 0) and space metric 
hik = diag(L^, i?^, i?^ sin^ 6') given foliation of the manifold on space and 
time. The scalar curvature density has the form 

N^/h{^^^R+{iTtKY -TiK^))- (5 3) 

-2 (yhK)^oj + 2 (VhKN' - Vkh'^j^ 

where ^^^R and K^^ are the scalar curvature of a space metric hij and exterior 
curvature tensor of a surface t = const. Substituting expression (|5.1| ) for the 
metric into ( |5.3| ) we obtain the expression for internal and external curvatures 
of the surface t = const in the form 

and 

irj = diag(ir;,ir,^KJ) 

= j^{l- L'N^ - L(iVO') , (5.5) 



Kl = K; = ^i^R-R'N^ 

Here dot and prime denote differentiation in t and r respectively. Contribu- 
tions to the gravitational action from the terms containing total derivatives 
in ( p.3|) give rise to the surface terms which cancel each other at the common 
boundaries of regions I, II and II, III. So we are left with the surface terms 
at infinity. We will turn to them later. The essential part of the action for 
gravitational field is just the ADM part of the action ( ^.2|) with Lagrangian 



V = TTTT^NLR^ ( (3)^ - (M)' - Tri^^) (5.6) 

Contribution to the action from the integral over the region III in the limit 
e — > is only due to the term containing second derivative of i?, namely 

/ 1 ^,,«M3,«_/^_/MM (5.7) 

J levrG J GL J GL 



III II 
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We will denote by hats variables on E and by [A] = lini^^o {A.{r + e) — A{r — e)) 
a jump of variable A{r) on the shell surface. Substituting the expression ( |5.1|) 
into the shell part of the action we have: 



Ssheii = -M J ^m-L^ (^Nr + f ) 'dt (5.8) 



The explicit form of the action ( |5.2| ) with metric ( |5.1| ) becomes 



The canonical formalism for this action can be described in the following 
way. Momenta conjugate to corresponding dynamical variables are 

= 6S/6N = 0; 

P^r = SS/5N^ = 

Pl = 6S /6L = -^ - 

Pr = 5S /5R = -^ I^R'N' ((^^'■)' - ^] 



p^ = 5S/5R = 
p^ = 6S/6L = 

TT = Ob / or = . ^ ' 
The action fB.OI) rewritten in the Hamiltonian form becomes 



(5.10) 



S= J (PlL + PrR -NH - N'-Rr) drdt + / 7rf- 

N (^R[R']/{GL) + ^f^^^^'^^yiPj - (5.11) 
[-L [Pl] - tt) dt 

with 

^ - ^\2W~^~)^'0\-^-^^\nr) ) (5.12) 
Hr = PnR'-LPi 
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where N, N^, N and are Lagrange multipliers in the Hamiltonian formal- 
ism. The system of constraints contain two surface constraints in addition 
to usual Hamiltonian and momentum constraints of the ADM formalism. 
ADM constraints: 

' H = 



(5.13) 



^ Hr = 

Shell constraints: 

■ Hr = n + L[PL]=0 
R[R' 



(5.14) 







Karel Kuchar proposed some specific canonical transformation of the vari- 
ables {R, Pr, L, Pi) to new canonical set (-R, Pr, M, Pm) in which Hamilto- 
nian and momentum constraints given by ( |5.12|) are equivalent to the very 
simple set of constraints : 

Pr = ^ 

(5.15) 

M' = 

The idea is to use the Schwarzschild anzatz for the space-time metric instead 
of the metric (|5J 



ds^ = -F(R, m)dT'^ + — ^ — -dR^ + R^idO'^ + sin^ Odcj)'^) (5.16) 

F(it, m) 

where T, R and m are some functions of (r, t) and F{R, m) = 1 — 2Gm/ R, 
and m, in general, is a function of r, m = m{r). Equating the two forms 
of the metric ( ^.11) and ( |5.16D we obtain the transformation between the two 
sets of dynamical variables. The explicit form of the transformation is 



CRT 

F 



-FPl 



p _ RFPjn {RY_ pp2 

- G~\/^ (5.17) 

R = R 

P - P , Pm , FP^ , {RFPJ'RR' -RFPmiRR'y 

- ^R+2n + ^+ GRF{{R!f / F-FPl) 
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where P„ 



-T'. The Liouville form 







PrR + PlL 



(5.18) 



can be expressed in the new variables as follows: 



e 



+ 



d 



Pmrh + PrR + 

1 



d 



dr \2G 



-RRln 



dt 

RR! + LPlG 



LPl + —RR' In 



RR' - LPlG 



RR' + LPlG 



rr'-lplG '• ^^-^^^ 

When there is no shell the total derivatives in ( p.l9| ) give rise to some surface 
terms at infinities. As was shown by K.Kuchar the appropriate falloff condi- 
tions at infinities make the last surface term in the Eqn.( prT9| ) zero. Then it 
follows from ( ^.191) that {R, Pr, m, Pm) form the canonical set of variables, 
and Eqn.( |5.l7| ) describes a canonical transformation between [R, Pr, L, Pl) 
and {R, Pr, m, P^) . 

In the presence of the thin shell the situation is different. Now surface 
the terms should not be neglected. Let us do the transformation ( ^.17| ) in 
regions I and II of our space-time. The Liouville form of our Hamiltonian 
system ( |5.11| ) has the form 



e 



PrR + PlL 



nr. 



(5.20) 



i+ii 



After integration the total derivatives in (|5.19|) give some contribution to the 
Liouville form on S: 



e 



PrR + P^m + 



i+ii 



LPL + —RR'\n 



RR' - LPlG 



RR' + LPlG 



rdt 



2G 



RRln 



RR' + LPlG 



RR' - LPlG 



+ 



Tcr 



P^m + PrR + pr + 



PrR 



(5.21) 



I+II 

where we denoted 



P 



n + L[PL] 



P, 



R 



± 



2n 



Pin 



RR' - GLPl 
RR' + GLPl 



(5.22) 
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and made use of the identity 



k = jR{t, m = m, r) + r)f (t)) |,^,(,) (5.23) 

The sign in front of logarithm in the definition of the momentum depends 
on whether we intersect the shell from "in" to "out" when going along the 
time coordinate curve from past to future ("+"-sign) or the other way around 
( "-"-sign). We see that this canonical transformation involves all the set of co- 
ordinates in the phase space 11 = {(-R(r, t), Pr^v, t), L{r, t), PL^r, t),f(t), 7r(t))} 
according to the formulae (^.17|) and ( ^.22| ). Moreover it introduces addi- 
tional pair of canonically conjugate variables {R, P^) on the shell. In both 
inner and outer regions I and II constraints are simplified due to the canon- 
ical transformation as it was in the absence of the shell (|5.15|) . The surface 
momentum constraint Hr = ( |5.12| ) takes the form 



p = (5.24) 

To go further we need to know the Hamiltonian constraint on the shell 
in terms of these new canonical variables. To do this, let us consider more 
carefully Eqn.( pT^3[ ) which is the full time derivative of the shell radius. Using 

the definition Pl = --^ (ji- we get: 

C N P 

R = -——^ + R'{f + N'-) (5.25) 
R 



Remembering that 

TT 



ML^ {N"^ + r) 
m - L2 (A^*- + rf 



we can find 

L [N' + r) 



Eqn.(|;2|) now reads 

R GPr R! TT 



N R L 



(5.26) 
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It is now easy to see that the jump of R across the shell is a linear combination 
of constraints 

GN ( „ H'' 



R 



R 



XH-—] (5.27) 



where 

vr 

X = ■ 



Now, the definition of the momentum R can be rewritten as follows 

(tl^i!il^v^,^±^i±4 (5,28) 

V / out 

where a = ^ and y = G^. 

The next step is to find the relation between a and y. From the definitions 
of Pl and R we have 



R!^ ^ , fay 
^ \X X 



. Solving this for y we get 



y = ^^V^'- - ^) - (5.29) 



Substituting this expression back into Eqn. (|5.28| ) we obtain after simple 
algebraic transformations 



P = ^^^v^g^ S. (5.30) 

Z - (ToutVz'^ + Font 

where = (5.31) 

1 - X 

Making use of the momentum constraint we can wright the jump of y as 
follows 



(JoutXV'^'^ + (1 -X^)Fout = OinX^Op' + (1 -X^)Fin - X^T^-X^ 

(5.32) 
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Note that this is just our starting equation ( P^ ) if we choose the proper time 
(substituting for z^.) Squaring this equation we get 



^in^J~^ ~l~ Pin ~ 
^out Fout 



RjF] , MG 
2MG 2R 



2MG 



MG 
2R 



(5.33) 



(5.34) 



The sign for z is opposite to the sign for R. Finally we get for the shell 
constraint 

P _ ^ 2MG 



MG 

2R_ 

y I _L MG_ 

^ 2MG ^ 2R 







(5.35) 



It can be easily shown that the above constraint is equivalent to the 
Hamiltonian constraint ( p.30[ ) . In what follows we will consider the following 



squared version of the Hamiltonian constraint (|5.3CI|) as the suitable classical 
counterpart for the quantum constraint for the wave function \1/ 



C — Font + Fi 



exp 



R 



R 



exp 



R 



i?2 

(5.36) 

The Hamiltonian constraint (|5.3CI|) was derived under the assumption that 
both Fin and Fout are positive . It is possible, of course to derive analogous 
constraints in T±-regions, where F < 0. But, instead, we make the following 
substitution 

^ F^/^ (5.37) 

and consider this function as a function of complex variable. Then the point 
of the horizon F = becomes a branching point , and we need the rules of 
the bypass. We assume the following 



pi/2 ^ |_p| e*'^' 

6 = in i?+-region 

- 7i/2 in T_ -region 

) = TT in i?_-region 

-n/2 in T+-region 



(5.38) 
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for the black hole case , and 



(p = 71 in i?_|_-region 

(f) — — 7r/2 in TL-region 

(f) — in i?_-region 

= 7r/2 in T+-region 



(5.39) 



for the wormhole case. The reason for such analytical continuation is that we 
are able to get the single equation on the wave function ^ which covers all four 
patches of the complete Penrose diagram for the Schwarzschild spacetime. 
The Carter- Penrose diagram for a = ±1 (balck hole and worm hole case, 
correspondingly) are shown in Fig. 10. 



R=0 




R=0 



R=0 




R=0 



Figure 10: 

It should be stress here that we consider the thin selfgravitating shell 
together with the space-time which it "creates" . Thus for different type of 
shell's motion we have to choose different forms of Hamiltonian constraints. 
For the shell escaping to infinity we choose the following constraint (see Fig. 8 
for the corresponding Carter-Penrose diagram in the case Fj„ = 1) 

V inV out 

while for the falling shell that start from infinity (Carter-Penrose diagram 
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is shown in Fig. 9) we have 



y/ fin's/ -P^out 

The bound motion is symmetric to the change of the momentum sign, so 
for the Hamiltonian constraint we choose the sum of the constraints (|5.40|) 
and (|5.41|) , namely 

e§^+e-^=a °"*tl!"^ (5.42) 

We now turn to the Dirac constraint quantization procedure . 

It is convenient to make a canonical transformation from {R, P^) to 

^ s 



i?2 



Rn 



Ps = Rl^ 



p° (5.43) 

LA. 
2R 



where Ro is some value radius of the shell. Dimensionless variable S is the 
surface area of the shell measured in the units of horizon area of the shell of 
mass M. 

The phase space of our model consist of coordinates {R{r), Pff{r),m{r), 
Pmij"), S, Ps, f,pr) f £ (— C)0, f — e) |J(f + e, oo). Then the wave function in 
coordinate representation depends on configuration space coordinates: 

^ = ■^{R{r),m{r), S, r) (5.44) 

and all the momenta become operators of the form 

pR{r) = -I 5/ 6R{r) P^r) = -i 5/ 5m{r) 

= -I d/ dS pr = ~id/df ^ ' ' 

ADM and shell constraints ( |5.13| ) and ( |5.14| ) become operator equations on 
\E'. The set of ADM constraints is equivalent to the set of constraints (|5.15|) 
in Kuchar variables which could be easily solved on quantum level. Indeed, 
in the regions I and II the equations 

d^/dR{r) = 

M'(r)^ = ^ ' 
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express the fact that wave function does not depend on R{r) and the depen- 
dence on M (r) is reduced in each region I and II to \1/ = 6{M — M±) where 
M± defined in the regions I (-) and II (+) do not depend on r. They equal to 
Schwarzschild masses in the inner and outer regions Mm and Mout in ( |4.19|) . 
The set of shell constraints ( ^.141) impose further restrictions on First 



of them takes the form 

5^/9f = (5.47) 

in new variables according to ( ^.24|) . So the only nontrivial equation is the 
shell constraint ( |5.36D (or ( p.30[ ( [OB| ) , ( [4. 191 ) which are classically equivalent 



to the Eqn.(|;36D) 

C (m+ ,m.,S, -ihd/dS^ =0 (5.48) 
M/ = vl/(m+,m_, S") 

The operator C contains the exponent of the of the momentum Pg. This 
exponent becomes an operator of finite displacement when Ps becomes dif- 
ferential operator: 



GPH £S ^ g 

e~ = e'^ovl/ = e *^ssv]/ = \[/(m+,m_,S' - (i) (5.49) 

A _ 2G 

For the sake of simplicity we consider here only the case of bound motion 
of the single shell having the bare mass M. It means that mj„ = and 
we should be use the symmetric form Hamiltonian constraint. In this case 
it is naturally to choose for Rq the Schwarzschild radius Rq = Rg = 2Gm 
(m = niout)- For ( we have now ( = |(^)^,and the Hamiltonian constraint 
can be written in the form 

2 1 m2 

vl/(5 + <) + vl/(5 - <) = ^' ^(S). (5.50) 

We have mentioned already that the classically equivalent constraints give 
inequivalent quantum theories. This is well known fact. We suggest that the 
criterion to choose the correct quantum theory is the behavior of the wave 
functions in the quasi classical regime. In our case this means the large black 
holes limit. Indeed , the parameter ( = |(-^)^ becomes small for large 
masses , and the expansion with respect to this parameter is equivalent to 
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the expansion in Planckian constant h {nipi = ^Jh/Gc). In the next Section 
we will consider this quasiclassical limit and show that our choice for the 
quantum constraint is a good one (at least in the case of one thin shell). 
At the and of this Section we would like to make an important remark. 
Our quantum equation ^.50| (which is just a Schroedinger equation ) is the 
equation in finite differences rather than differential equation, and the shift in 
argument is along an imaginary axis. In the case of differential equation we 
we require the solution to be differentiable sufficiently many times . Similarly, 
we have to demand the solutions of our finite differences equation ^.50| to be 
analytical functions. This condition is very restrictive but unavoidable. Our 



previous experience (see |^) shows that it is the analyticity of the wave 
functions and not the boundary conditions that lead to the existence of the 
discrete mass (energy) spectrum for bound states. 

The finite difference equation ( |5.5CI| ) becomes an ordinary differential 
equation in quasiclassical limit which is the same as the limit of large ( 
m ^ nipi ) black holes. Indeed the parameter of finite displacement of the 
argument of \I' in ( p.50| ) C, becomes small and we could cut the Tailor expan- 
sion 

^{S + CO = ^(5) + <^'(5) - ^^"{S) + ... (5.51) 
at the second order. Thus, 



(5.52) 

This leads to ordinary differential equations of second order, which are dif- 
ferent in _R_, T+ and T_ regions due to the different values of the phases 
in Eqn.( p750| ). The interesting for us singular points of these differential 
equations are 

S = oo and S = 1. (5.53) 

In the quasiclassical limit our requirement of the analyticity of the solutions 
to the exact equation (|5.50|) transforms into the requirement that the branch- 



ing points of the leading terms in the solutions to the approximate equations 
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should be of the same kind . Thus, we need to keep only those terms in the 
corresponding equations that give us these leading terms. Below we consider 
the black hole case only. The results are easily translated to the wormhole 
case. 

The singular point 5* = oo in the region lies in a classically forbidden 
region as far as we restrict ourselves with bound motions of the shell only. In 
order to analyze the behavior of \1' in this region we should take (|5.5CI| ) with 
= and expand all the quantities in terms of where s = (1 + yY The 
result is 

1 1 / 1 \ 

- + 72 1 - — + ^(2 - — ) P = (5.54) 
y Q \ m 2,y m J 

The leading term of the solution is 

m 



\ij ^ y2 AfiC exp{-iJ,y), (5.55) 

1 [a? " 
C V 

For another singular point in i?+ region , that is for 5 ^ 1 + we have 



^^fr^l-"{l- ) ^3/2 (5.57) 



with leading term 



3C^ V 

y = v^, s > C, y > C, C < 1 

Comparing the types of the branching points at s — oo and s ^ 1 + 
we can conclude that 

^ 2 

^ — n, n = integer (5.58) 
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This is the first quantization condition. We will not consider here the worm- 
hole case . Note only that , as can be shown , the positive values of quantum 
number n correspond to black holes while negative n correspond to worm- 
holes. 

We do not consider here separately the asymptotics in i?_ region near 
the horizon (s — > 1 -|- 0) because it differs from the corresponding solution in 
i?+-region only by the sign in front of the second term. 

Let us now turn to the asymptotics of the solutions in i?_-region for 
s — > cxo. Due to the minus sign in front of F^/^ the equation for the wave 
function in a i?_-region is quite different from that in a i?_|_-region 

1 1 / \ 

- - 72 1%' + 1 - — ) U = (5.59) 

y C \ m J 

The leading term of the asymptotic is now the following 

2 

m 

,2 



8C exp-(-y^) (5.60) 

Note that the falloff in the i?_-region is much faster than it is in the 
region. This is a quite reasonable result because it means that the quantum 
shell in the black hole case can penetrate into the it!_-region (which is com- 
pletely forbidden for the classical motion) but the probabihty of such an 
event is negligible small. 

And , again, comparing the types of the branching points at s ^ 1 -|- 
and s — > cxD in the i?_-region we get 



— -1 1 

^ — = - + p, p — positive integer (5.61) 

8C 2 

The appearance of the second quantum number is rather surprising and its 
origin is in nontrivial casual structure of the Schwarschild space-time. The 
classical shell in the (e.g.) black hole case cannot move in the i?_-region, and 
simply absent on the complete Carter-Penrose diagram. But the quantum 
shell with parameters corresponding to the black hole case, "feel" both 
and i?_-regions, this why the second quantum number plays its role. It is 
instructive to compare our result with the famous Coulomb problem for the 
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hydrogen atom. The parameters of the Coulomb problem are the proton's 
and lectron's electric charge and all of them being quantized from very begin- 
ing. Therefore, to obtain a discrete spectrum for bound states it is enough to 
have only one quantum condition. But, in our case there are two continous 
parameters, the total mass m and bare mass M. And now, in order to obtain 
a discrete mass (energy) spectrum, we should have two quantum conditions. 
And just what we really got! 

Now, we would like to consider the behavior of the solutions in the vicin- 
ity of the horizon (sub-Planckian deviation ), where \y\ 3> C ("5 ~ !)• To 
be specific we will be interested in the solutions i?+ and T_ regions. The 
expansion ( 5.52 ) is no more valid for the function F^/^(s ± iQ but it is still 



valid for the wave \E'. Keeping the leading terms only we have now 
2 , , / 4 , M\ 2 



vl>,,(.) - -^,{s) + { -:^(1 - -^) - ) vi/(s) = (5.62) 



with the solution 



The coefficient a equals to 1 in the i?_|_-region and to imaginary unit i in the 
T_ -region. 

In the i?+ region 





(5.64) 



and we have superposition of two waves (ingoing and out outgoing) with 
relatively equal amplitudes, 
in the T_ -region 



^){l + z)= (5.65) 




4m^' 
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The existence of two waves in the T_-region reflects the quantum trembling 
of the horizon . But the outgoing wave is enormously damped relative to the 
ingoing wave (of course , in the T_- region the situation is exactly inverse 
one). It is this damping that cause (in a quasiclassical regime) existence of the 
single ingoing wave in the T_-region at the distances larger than Planckian. 

In the end of this Subsection let us discuss briefly some features of the 
obtained discrete mass spectrum (|5.58| ) and (|5.61|) 



Let us discuss some properties of the spectrum that arises from these 
conditions. 

1. For larger values of quantum number n (^ — 1 << 0) one can 
easily derive nonrelativistic Rydberg formula for Kepler's problem, Enonrei = 
M-m = -^. 

2. The role of turning point po is now played by the integer n. Thus, 
keeping n constant and calculating 7 = |^|n one can distinguish between 
a black hole case (7 > 0) and a wormhole case (7 < 0). It appears that 

> for n > riQ, negative or zero, and 

l^ol = E[V2\I - 2%l + 2j9)] (5.66) 

where E[- ■ ■] means the "integer part of It is interesting, that for n > 
the asyptotical behaviour of the wave function is exactly the same as in the 
Coulomb problem, but for the negative values of no there is no such analogy. 
3. There exists a minimal possible value for a black hole mass. This occurs 
if p = no = 0, 

mjnin = V2mpi (5.67) 

4. It is interesting to study the limiting case when M — >■ 00 and total mass 
m is kept finite. It is correspond to the limit of eternal black hole. From our 
spectrum it easy to deduce that 

m = lim ( ) mpi, (5.68) 

n,p^oo \^ ^ J 

and, as result, we obtain that mass of the black hole may be arbitrary (but 
positive). 

5. The our spectrum is not universal in the sense that corresponding wave 
function form two-parameter family \l/„,p. But for quantum Schwarzschild 
black hole we expect a one-parameter family of wave functions. Quantum 
black holes should have no hairs, otherwise there will be no smooth limit to 
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the classical black holes. All this means that our spectrum is not a quantum 
black hole spectrum, and our shell does not collapse (like an electron in 
hydrogen atom). Physically, it is quite understandable, because the radiation 
is yet included into consideration. But before doing this we would like to 
make one important note. 

Till now we considered only the case of zero mass inside the shell. The 
generalization to the general case is not difficult, and result is the following 



2(Am)2 - 



2mpin, 



v/M2 - (Am)2 

(5.69) 

- (Am)2 = 2m^;(l + 2p), 

Where M is the bare mass of the shell, m = niout is the total mass of the 
system, and Am = niQut — friin. Moreover this spectrum is, in fact, the ex- 
act spectrum of our finite difference equation. The reason why the formulas 
obtained for large black holes appeared the exact ones is pure mathemati- 
cal and we will not discuss it here. But one important feature of the exact 
spectrum (|5.69|) should be mentioned. We know already that all linear ho- 
mogeneous equation in finite differences obey the following feature, namely, 
any its solutions multiplied by the periodic function which period equals to 
the argument shift of the equation, gives again the solution. This leads, in 
general, to infinite degeneracy of the spectrum. But our spectrum ( |5.69| ) 



can be made nondegenerate. The matter is that as the argument of our 
equation in finite differences ([5. 50]) we have square of the radius and leading 



term of the asyptoticalsolution at the right infinity (in i?+-region), in general, 

2 

has a form \E' ~ e"^ , but there is one exception. The exceptional solution 
has an asymptotic \Ef ~ e^^, and just the same asymptoticly behavior ex- 
hibit the wave function of the nonrelativistic Coulomb problem. Thus if we 
demand the existence of the "correct" non-relativistic limit, then "correct" 
wave function should have the leading asymptotical term ~ e^^', and such 
a requirement removes a degeneracy. 



5.1 Radiation and quantum collapse dynamics 

In order to investigate the quantum collapse, we should switch on the process 
of quantum radiation. For the modelling of such radiation we again make 
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use of the developed above theory of thin shells, but now shells are null and 
move with the speed of light. The corresponding classical equation is our 
nonsymmetric constraint ( p.40| ) in which the bare mass M is put zero. The 



quantum equation looks as follows p5 | 



1- 

_ _^ '^{min, mout, s) (5.70) 



here /x = rriin/ rriout, C = \^pi/''^ouf The existence of the second infinity on 
the other side of the Einstein-Rosen bridge leads to the following quantization 
condition (m = rriout) 



5m = mout - min = -2m + 2./m2 + kml^, (5.71) 



where k is an integer. It is interesting to note that if we put k = 1 (minimal 
radiating energy) and require 6m < m (not more than the total mass can be 
radiated away), then we obtain 

2 

m = mout > —j=mpi. (5.72) 
V5 



Thus, the black hole with the mass given by Eqn. ( ^.67|) is not radiating and. 



therefore, it can not be transformed into semi-closed world (wormhole-like 
case). 

The discrete spectrum of radiation ( ^.71j ) is universal in the sense that 
it does not depend on the structure and mass spectrum of the gravitating 
source. This means that the energies of radiating quanta do not coincide 
with level spacing of the source. The most natural way in resolving such 
a paradox is to suppose that quanta are created in pairs. One of them is 
radiated away, while another one goes inside. Thus, the quantum collapse 
can not proceed without radiating even in the case of spherical symmetry. 
This radiation is accompanying with creation of new shells inside the primary 
shell we started with. We see, that the internal structure of quantum black 
hole is formed during the very process of quantum collapse. And if at the 
beginning we had one shell and knew everything about it, then already after 
the first pulse of radiation we have more than one way of creating the inner 
quantum. So, initially the entropy of the system was zero, it starts to grow 
during the quantum collapse. If somehow such a process would stop we would 
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call the resulting object "a quantum black hole". The natural limit is the 
transition from black hole to the wormhole-like shell. The matter is that 
such a transition requires (at least in quasi-classical regime) insertion of an 
infinitely large volume, and the quasi-classical probability for this process is 
zero. 

Let us write down the spectrum of the shell with nonzero Schwarzschild 
mass, the total mass inside, mj„ 7^ 



2 (Am) 2 - M2 2m, 



- (Am)2 Am + m. 



n (5.73) 



- (Am)2 = 2(1 + 2p)m^ 



Here Am is the total mass of the shell, M is the bare mass, the total mass 
of the system equals m = irtout = + mj„. For the black hole case 
< AmAm, or 



^>i(/(^PH-l-^). (5.74) 

After switching on the process of radiation governed by Eqn. (|5.7T| ), the quan- 
tum collapse starts. Our computer simulations shows that evolves in the 
"correct" direction, e.g. it becomes nearer and nearer to the threshold ( |5.74| ) 
between the black hole case and wormhole case. The process stops exactly 
at n = 0! 

The point n = in the spectrum is very special. Only in such a state 
the shell does not "feel" not only the outer regions (what is natural for the 
spherically symmetric configuration) but it does not know anything about 
what is going on inside. It "feel" only itself. Such a situation reminds 
the classical (non-spherical) collapse. Finally when all the shells (both the 
primary one and newly produced) are in the corresponding states rii = 0, the 
system does not "remember" its own history. And this is a quantum black 
hole. The masses of all the shells obey the relation 

Ami = -^Mi. (5.75) 

The subsequent quantum Hawking's evaporation can produced only via some 
collective excitations and formation, e.g., of a long chain of microscopic semi- 
closed worlds. 
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5.2 Classical analog of quantum black hole 

Let us consider large (m ^ rripi) quantum black holes. The number of shells 
(both primary ones and created during collapse) is also very large, and one 
may hope to construct some classical continuous matter distribution that 
would mimic the properties of quantum black holes. First of all, we should 
translate the "no memory" state (n = for all the shells) into "classical lan- 
guage". To do this let us rewrite the Eqn.( p.7|) (energy constraint equation) 
for the shell, inside which there is some gravitating mass rriin, 




2Gm.;r, _ 2Gmn„t GM 



P 

and consider a turning point, {p = 0, p 




(5.76) 



Am = rriout - rriin = M i 1 — . (5.77) 



It is clear now that in order to make parameters of the shell ( Am and M) 
not depending on what is going on inside we have to put rriin = opo- 

Our quantum black hole is in a stationary state. Therefore, a classical 
matter distribution should be static. We will consider a static perfect fluid 
with energy density e and pressure p. A static spherically symmetric metric 
can be written as 

ds^ = e^de - e^dr^ - r^dO^ + sin^ Odip^) (5.78) 

where u and A are functions of the radial coordinate r only. The relevant 
Einstein's equations are (prime denotes differentiation in r) 

8nGe = -e\\--) + \, 

-SttGp = -e^{\ --) + \, (5.79) 

-S^Gp=--e>-W' + - + — -) 

The first of these equations can be integrated to yield 

g-^ = 1 (5.80) 
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where ^ 

m{r) = An er''^dr' (5.81) 



is the mass function, that must be identified with mj„. Thus, m(r) = ar, 
and 



e = e-^ = l- 2Ga. (5.82) 

We can also introduce a bare mass function 

M{r) =A7r eeh''^dr', (5.83) 
Jo 

and from Eqn.( ^.8^ ) we get 

The remaining two equations can now be solved for p{r) and e'^. The solution 
for p(r) that has the correct non-relativistic limit is 

1 

47rr2' " G' 



p(r) = ^— ^, 6 = -(1 - 3Ga - Vl - 2GaVl - 4Ga), (5.85) 



and for e"^ we have 

e'' = Cr^^^^. (5.86) 

The constant of integration G can be found from matching of the interior and 
exterior metrics at some boundary r = vq. Let us suppose that r > vq the 
space-time is empty, so the interior should be matched to the Schwarzschild 
metric. Of course, to compensate the jump in pressure {Ap = p{ro) = po) 
we must introduce some surface tension E. From matching conditions (see, 
e.g. [Ill) it follow that 

_2G-Ji+b_ 
2G, 



G = (1 - 2Ga)ro 



e'' = (l-2Ga)(-)'^T^, (5.87) 
^ 2Ap b 



3 



ro 27rr^ 

We would like to stress that the pressure p in our classical model is not real 
but only effective because it was introduce in order to mimic the quantum 



66 



stationary states. We see, that the coefficient b in Eqn.( p.85| ) becomes a 
complex number if a > 1/AG. Hence, we must require a < 1/4:G, and 
in the hmiting point we have the stiffest possible equation of state e = p 
It means also that hypothetical quantum collective excitations (phonons) 
would propagate with the speed of light and could be considered as massless 
quasiparticles. It is remarkable that in the limiting point we have m(r) = 
M(r) / - the same relation as for the total and bare masses in the "no 
memory" state n = 0! The total mass mo = m{rQ) and the radius tq in this 
case are related mo = 4G'ro - twice the horizon size. 

Calculations of Riemann curvature tensor R^i^ and Ricci tensor Rik show 
that if p < e {a ^ b) there is a real singularity at r = 0. But, surprisingly 
enough, both Riemann and Ricci tensors have finite limits at r — 0, if e = p 
(a=b=l/4G). Therefore we are allowed to introduce the so-called topolog- 
ical temperature in the same way as for classical black holes. The recipe 
is the following. One should transform the space-time metric by the Wick 
rotation to the Euclidean form and smooth out the canonical singularity by 
the appropriate choice of the period for the imaginary time coordinate. The 
imaginary time coordinate is considered proportional to some angle coordi- 
nate. In our case the point r = is already the coordinate singularity. The 
azimuthal angle has the period equal to vr. Thus, all other angles should be 
periodical with the period vr. The topological temperature is just the inverse 
of this period. 

The easy exercise shows, that the temperature 

T=7^ = = Tbh (5.88) 

exactly the same as the Hawking's temperature Tbh 0! The very possibility 
of introducing a temperature provides us with the one-parameter family of 
models with universal distributions of energy density and pressure 

the parameter being the total mass uiq or the size tq = 4Gmo. It should be 
noted that the two-dimensional part of the metric (|5.78|) in the limiting case 



a = 6 is nothing more but the Rindler's metric, and the null surface r = 
serves as the event horizon. The Carter-Penrose diagrams for this case is 
shown in Fig. 11. 
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Figure 11 
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We can now develop some thermodynamics for our model. First of all we 
should distinguish between global and local thermodynamic quantities. The 
global quantities are those measured by a distant observer. He measures the 
total mass of the system mo and the black temperature Tbh = T^o and does 
not know anything more. Let us assume that this observer is rather educated 
in order to recognize he is dealing with a black hole and to write the main 
thermodynamic relation 

dm = TdS. (5.90) 

In this way he ascribes to a black hole some amount of entropy, namely, the 
Hawking-Bekenstein value ^ 

S = i(i^I^ = 4.G™? = M^)= (5.91) 

4 Mpi 

The local observer who measure distribution of energy, pressure and local 
temperature is also rather educated and writes quite a different thermody- 
namic relation 

e{r) = T(r)s(r) — p{r) — fi{r)n{r). (5.92) 

Here e{r) and p{r) are energy density and pressure, T(r) is the local temper- 
ature distribution, s(r) is the entropy density, /i(r) is the chemical potential, 
and n(r) is the number density of some (quasi)" particles" . For the energy 
density and pressure the local observer gets, of course, the relation ( ^.89|) , 
and for the temperature - the following distribution 

Tir) = (5.93) 
V /vrr 

which is compatible with the law T(r)e5 = const and the boundary condition 
Too = Tbh- Such a distribution is remarkable in that if some outer layer 
of our perfect fluid would be removed, the inner layers would remain in 
thermodynamic equilibrium. And what about the entropy density? Surely, 
the local observer is unable to measure it directly but he can receive some 
information concerning the total entropy from the distant observer. This 
information and the measured temperature distribution ( |5.93| ) allows him to 
deduce that 

s(r) = — ^ — (5.94) 

and 
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It is interesting to note that in the main thermodynamic equation the 
contribution from the pressure is compensated exactly by the contribution 
from the temperature and entropy. It is noteworthy to remind that the 
pressure in our classical analog model is of quantum mechanical origin as 
well as the black hole temperature. And what is left actually is the dust 
matter we started from in our quantum model, namely, 

We may suggest now that the quantum black hole is the ensemble of some 
collective excitations, the black hole phonons, and n{r) is just the number 
density of such phonons. 

Knowing equation of state, e = p, we are able to construct all the ther- 
modynamical potentials for our system. As an example we show here how to 
calculate the energy as a function of the entropy S, and the number particles 
N. By the first law of thermodynamics 

dE = TdS - pdV + i^dN (5.97) 

where T = |^ | ^ ^ is a temperature p — |^ | ^ ^ is a pressure, and n — §w\sv 
is a chemical potential. The energy is additive with respect to the particle 
number N, hence, E = Nf{x,y) where x = and y = y. Since e — ^ — 
yf{x, y) and p — y^^ from the equation of state we obtain 

/ = a{x) — na{x) 
e — p — n'^a{x) 

Further, 

T = na'{x) 
II — n{2a — xa') 

But, in any static gravitational field T — Tq/^/q^ and /x = fio/y/goo, so 
II — ^qT, where 70-some numerical factor. Thus, 

2a — xa' = 7oa', 
a{x) = Co(7o + xf 

70 



(5.96) 



where Cq is a constant of integration. It is easy to see that p/T — 1/ACq. In 
our specific model p/T^ = tt/SG, so Cq = 20/71. Moreover, because of the 
relation e = p = Ts = /xn we know that free energry F = E — TS is zero. 
From this we have for the entropy 

S = 7oiV (5.98) 

The black hole entropy equals one fourth of the dimensionless horizon area, 
and form this we recover the famous Bekenstein-Mukhanov mass spectrum 

J^^/Nrupi (5.99) 
V 47r 

We can even calculate the remaining unknown coefficient 70 using the 
phonon model. Indeed, since the free energy F = and from the well know 
relation 

F = -rin^e-^ (5.100) 

n 

we deduce that the partition function 

Z = J]e-^ = l (5.101) 

n 

Assuming now, that our gravitational phonon are linear oscillators, we have 

e 2T 

Z^- ^ = 1 (5.102) 

1 — e T 

and 

\/5 - 1 , , 

e"^ = — (5.103) 

The energy level is nothing more but bare mass spectrum. For the black hole 
mas spectrum we have (n ^ 1) 

m = V2Tn\J^ = ^i^^ (5.104) 

Identifying the level number n with the "particle" number N in our thermo- 
dynamical model, we obtain finally 

7o = In — ^ 0.5 

and for the black hole mass spectrum 

m « Q.2\fNmpi (5.105) 
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6 Conclusion 



And again, the question which aheady have been asked in Introduction. 
Why should we study black holes, for what? I hope, at least something 
became clear while reading this paper. The main feature of the black holes 
is their universality. It seems it is this that should distinct quantum black 
holes from other quantum objects, irrespective of what kind of theory the 
future quantum gravity will be - fundamental, effective, induced or some 
other. In the absence of elaborated quantum theory we have to construct 
models. It is desirable that such models were exactly solvable, because the 
exact solutions serve as a background of our physical intuition, which, in 
turn, allow us to separate the model depending features of the object under 
consideration (in our case - quantum black holes) from the true fundamental 
ones. And, in general, the author believes strongly, that in the future the 
specific features of quantum black hole will play the role of some kind of 
selection rule for selecting the "correct" quantum theories. It is like the 
"great Fermat theorem" in mathematics, where numerous attempt to prove 
it resulted in creation of new fields in the number theory. So,the investigation 
of quantum black holes is not simply the exciting intellectual game. 
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